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Abstract. We investigate, both analytically and by computer simulation, the 
kinetics of a microscopic model of hard rods adsorbing on a linear substrate. For 
a small, but finite desorption rate, the system reaches the equilibrium state very 
slowly, and the long-time kinetics display three successive regimes: an algebraic one 
where the density varies as 1/t, a logarithmic one where the density varies as l/£n(t), 
followed by a terminal exponential approach. A mean-field approach fails to predict 
the relaxation rate associated with the latter. We show that the correct answer 
can only be provided by using a systematic description based on a gap-distribution 
approach. 



In many situations, thermal energies are significantly smaller than the energy 
needed for hard particles to diffuse, and under applied external forces the 
system evolves towards non-equilibrium configurations, metastable configurations, 
or eventually reaches a stable equilibrium state after a very slow process. For 
example, the adsorption of some proteins and colloidal particles on solid surfaces Q 
involves particle-surface energies that are so strong that the process is characterized 
by extremely small desorption and surface diffusion on the experimental time scale. In 
granular materials, particles are trapped in a metastable configuration, unless external 
energy is brought to the system. Recent experiments^], [| measured the densification of 
a vibrated granular material. A column containing monodisperse spherical beads was 
tapped periodically with a given intensity and the powder evolved slowly, essentially 
as the inverse of the logarithm of the number of taps, from a loosely packed state to 
a denser steady state whose density depends on the tapping strength. In all these 
cases, geometric exclusion effects dictate the kinetics of densification, i.e., addition of 
new particles is exponentially limited by the inverse of the free volume |B |j. These 
effects are accounted for in a simple adsorption-desorption model || g, g. Its one- 
dimensional version, also known as the parking lot model, was shown to display a 
l/ln(£) approach to the final state for vanishing desorption rate|5|, |(|. More recently, 
Ben Nairn et proposed an approximate solution of the model for a small but 

non-zero desorption rate. We discuss here the solution and show that the mean-field 
treatment (or adiabatic approximation) used in the late stage of the process is not 
valid and greatly underestimates the characteristic time of relaxation. We develop 
a more consistent approximation which shows good agreement with the numerical 
simulations. 

The general definition of the adsorption-desorption model is as follows. Attempts 
are made to add objects to a space of arbitrary dimension at randomly selected 
positions with a constant rate k + . If the trial position does not result in an overlap 
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with a previously placed object, the new object is accepted. In addition, all objects in 
the system are subject to removal (desorption) at random with a constant rate In 
the parking lot model, the substrate is a line and the objects are hard rods. This 1-d 
model has been solved in some limiting cases. When fc_ = 0, the adsorption is totally 
irreversible and the process corresponds to a 1-d Random Sequential Adsorption 
(RS A) for which the kinetics are known exactly || . Due to the absence of relaxation 
mechanism, this process evolves towards a non-equilibrium state and the long-time 
kinetics are given by an algebraic scaling law, poo — pit) ~ 1/t, with ~ 0.747... 
When k + = 0, one recovers a desorption process for which analytical solutions are 
also available^). In the limit fc_ — > 0, accurate descriptions have been obtained^, 
In this case, the process cleanly divides into two sub-processes. The initial phase 
consists of an irreversible adsorption and it is followed by an infinite sequence of 
desorption-adsorption events in which a rod detaches from the surface and the gap 
that is created is immediately filled by one or two new rods. The latter possibility 
causes the system to evolve continuously to the close-packed state with p = 1, as || 
1 — pit) ~ 1/ ln(f ln(t)) where t now represents a rescaled time. For the general case, 
where both fc+ and fc_ are non zero, no complete solution is available. 

The properties of the parking lot model depend only on the ratio K = k + /k-. 
With an appropriate rescaled time, the kinetic equation describing the evolution of 
the density of adsorbed particles is given by 

where $, the insertion probability at time t (or density p), is the fraction of the 
substrate that is available for the insertion of a new particle. Thus, large values of K 
correspond to small desorption rates. The presence of a relaxation mechanism, even 
infinitesimally small, implies that the system eventually reaches a steady state that 
corresponds to an equilibrium configuration of hard particles with p cq — K$ cq (p eq ), 
where p eq denotes the equilibrium density. At equilibrium, the insertion probability 
is given exactly by 

$eq(p) = (l-p)e*p(-/V(l-p)). (2) 

Inserting Eq. (Q) in Eq. ([!]) leads to the following expression for the equilibrium density: 

P ^-\ + L w {K) (6) 

where L w {x) (the Lambert- W function) is the solution of x = ye y . In the limit of 
small K the isotherm takes the Langmuir form (p eq ~ K/ (1 + K)) while for large K, 
p eq ~ 1 — 1/ ln(AT). At small values of K, equilibrium is rapidly obtained, but at large 
values the densification is extremely slow. In the following, we restrict ourselves to 
the latter case. In the initial stages of the process, desorption events are negligible 
compared to adsorption and the process displays an RSA-like behavior. Figure 0a 
shows simulation results (see below for details) for the density versus 1/t in the 
intermediate density region (0.4 < pit) < 0.75). The larger the value of K, the 
larger the scaling law region. For p(t) > 0.7, which implies K > 100, desorption 
can no longer be ignored, because it permits particle rearrangements on the line and, 
eventually, insertion of additional particles. The mechanism for densification is similar 
to that of the model with infinitely small desorption and the kinetics is dominated by 
the l/ln(i) term (Figure Up). For large but finite values of K : this regime persists 
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until the density is very close to the equilibrium value and the desorption term is 
comparable to the adsorption term. In the final regime, an exponential approach to 
equilibrium is observed (Figure ^|a). This terminal approach corresponds to the linear 
response regime. 

To determine the kinetics of the densification process, a knowledge of &(p) is 
required. The evolution equations can be expressed by means of the gap distribution 
functions G(h, t) which represents the density of voids of length h. The time derivative 
of G(h, t) is given by 



dG(h,t) 
dt 



-H{h- l)(7i- l)G{h,t) +2 / dh'G(h',t) 

Jh + l 

- ^G(h, t) + H %~^ £ 1 dh'G(h', h-l-ti, t), (4) 

where G(h,h',t) is the gap-gap distribution function and H{x) is the unit step 
function. The first two terms on the right-hand side of Eq. (|4|) correspond to the 
loss and gain terms due to adsorption while the remaining two arc due to desorption. 
A similar equation can be written for G(h, h! ,£), and it requires higher-order gap 
distribution functions. As expected, even though it is one-dimensional, this kind of 
process leads to an infinite hierarchy of equations involving an infinite set of multi- 
gap distribution functions. Only approximate solutions can be obtained, by means of 
closure ansatz that enables the truncation of the hierarchy. 

The insertion probability can be expressed in terms of the gap distribution function 

as 



= J dh(h-l)G(h,t), (5) 
and we have in addition the following sum rules: 



p(t) = dhG(h,t) = l- dhhG{h,t). (6) 
Jo Jo 

The exact expression for the equilibrium structure is known ||, |) ; the gap distribution 
function is 

P 2 _ ( P 



G eq (h, p) = exp - -*—h , (7) 

l-p \ 1 - p J 

and all higher-order distribution functions satisfy the factorization property, 

G cq (hi,h 2 , ■ ■ ■ ,h„,p) = G cq (hi 1 p)G cq (h 2 ,p) ■ ■ -G cq (h n ,p). (8) 

An adiabatic (mean-field) approximation assumes that, at any density p(t), the 
structure of the adsorbate acquires very rapidly an equilibrium form satisfying Eqs. (0) 
and (||). This leads to an expression for $ akin to Eq. (g) with p(t) in place of 
p e g. Therefore, expanding Eq. (|l|) to first order in density, 5 pit) = p(t) — p^, with 
Poo = Peq{K), one obtains 

j t 5p=-r MF (K)5p + 0(Sp 2 ) (9) 

with 

(l + L w {K)f 

1 Mf{K) = — , (10) 
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which implies an exponential approach to the equilibrium state with a relaxation time 
given by K / \n{K) 2 for large K . Using a very efficient algorithm for the adsorption- 
desorption processes JTcj ] , we have performed Monte-Carlo simulations for K ranging 
from 10 to 5000. Averages were taken over 5000 independent runs. Figure |^b shows 
the relaxation rate versus K: the full curve gives the mean- field prediction, Eq. (|l0|), 
and dots correspond to simulation results. It is evident that the mean-field analysis 
gives a poor estimate of the relaxation rate for large K . In contrast with adsorption 
models with fast diffusion []l2[, the mean-field approximation is not valid here; the 
desorption mechanism becomes so inefficient for large K that the system does not 
follow a quasi-static path even for the regression of small fluctuations to equilibrium. 
In Figure one observes that the insertion probability $>(t) displays a non-monotonic 
behavior at long times for K > 100, and approaches the equilibrium value from below. 
Hence, at long times, the time derivative of $ becomes positive. As the density 
remains an increasing function of time along the entire process, the density derivative 
of the insertion probability is positive and bounded by 1/K by virtue of Eq. ([!]). 
Consequently, this shows also that the relaxation rate is smaller than 1/K. (Let 
us recall that both the density and time derivatives of $ are always negative in the 
adiabatic approximation.) 

To obtain the leading term in the exponential approach towards equilibrium, when 
K is very large (but finite), we assume that, as for the steady state, \G(h,t)\ ~ 
exp(— IT/i), with n ~ \nK. As a consequence, if one defines p n (t) — J™ +1 G(h,t)dh 

and <& n {t) = f™ +1 (h — l)G(h,t)dh, then p n ~ $„_i ~ K~ n , so that if one looks for 
the dominant behavior in 1/K, it is sufficient to consider the first intervals in h. As in 
the adiabatic approximation, one can expand the gap densities in power of Sp(t) and 
keep only the linear term which gives rise to the exponential decay. More specifically, 
by introducing 

A(h\ dhxG(h,p) 
A{n) = poo 



(11) 

d P Poo 

where A(h) is a piecewise continuous function, Eq. (^) can be rewritten for < h < 1 
as 

2 -^ L A{h) = ^ / + °° dh'e- p -<- h '-»AQ. + h') (12) 
K K Jh 

where Poo = Poo/0- — Poo) is the equilibrium pressure for poo = Peq and 7 = — KT = 
-KSp(t)/Sp(t)\ Poo . FromEqs. (0) and one obtains 



pOO 

7 = 1 - P£, / dhhe- p ~ h A(l + h) (13) 
Jo 

whereas the sum rules in Eq. (|^) give 

/>oc />oo 

Poo / dhe- p °° h A{h) = -P^ / dhhe- p ™ h A(h) = 1. (14) 
Jo Jo 



When integrating the two sides of Eq. ([14]) between and 1, one obtains to first order 
in l/K 



(2- 7 ) = 2P£ j\ 
Jo 



dhhe- r °° h A(l + h) + 0(l/K). (15) 
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Combining Eq. ( |l3| ) with Eq. fll5| ) yields 7 = 0(1/ K). Hence, the relaxation rate is an 
order of magnitude smaller than predicted by the mean- field treatment. It is actually 
given by (recall T = j/K) 

r = { m - P ™ f dhhe ~ P °° kA ( 2 + h )) + 0(1/K 3 ). (16) 
T(3T(3 (17) 

In order to obtain explicitly the coefficient of the leading term in l/K 2 , one needs 
to solve Eq. (|j) for h > 1. This can be achieved by assuming that the factorization 
property for the two-gap distribution function is valid to a 0(1/ K). Doing this, and 
after some tedious algebra, we have been able to derive the expression for A(h), i.e., 
for the gap distribution function at linear order in Sp. (The result, which is too 
lengthly to be presented in this letter, is of course different from that of the mean- 
field approximation and it compares very well with the simulation datajl^].) It leads 
to the following dominant behavior for very large K: 

r - ^ (-) 

In Figure ||b the dashed curve corresponding to Eq. ([IJ) gives a good agreement with 
the simulation results, contrary to the mean-field predictions. 

We have also examined the time correlations of the density fluctuations Sp(t) 
around the equilibrium state. Specifically, we have evaluated 

c(t) = ^£ Sp(t')Sp(t + t')dt'. (19) 

We have verified that at long time, the correlation function decays exponentially with 
a relaxation rate which is the same as that obtained by Eq. (|l8|), in agreement with 
Onsager linear regression principle. Both the fluctuation-dissipation relation and 
the time translational invariance are satisfied. However, (interrupted) aging[[ll] is 
expected in the l/Zn(i)-regime. 

We thank R. Muralidhar who proposed the simulation algorithm. Travel support 
from NSF and CNRS is also gratefully acknowledged. 
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Figure 1. (a) Evolution of the adsorbed density versus 1/t for several values of 
K, (K = 100,500,5000 from bottom to top) and for a RSA process (K = oo, 
dashed curve), (b) Evolution of the adsorbed density to its equilibrium value. At 
intermediate times, the adsorbed density evolves as l/ln(t), while very close to 
equilibrium the density relaxes exponentially. 
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Figure 2. (a) Final exponential approach of the density to its equilibrium value, 
(b) Relaxation rate T for the approach to equilibrium. Upper curve: prediction from 
mean-field approximation, Eq. (|lo|). Intermediate curve: prediction from Eq. (|l8|). 
Open circles: from the numericalsimulations. 
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